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Abstract. The temperature dependence of bottomonium states at temperatures above and below Tc is presented, using non- 
relativistic dynamics for the bottom quark and full relativistic lattice QCD simulations for two light flavors on a highly 
anisotropic lattice. We find that the S-waves (Y and r/^) show little temperature dependence in this range while the P wave 
propagators show a crossover from the exponential decay characterizing the hadronic phase to a power-law behavior consistent 
with nearly-free dynamics at approximately twice the critical temperature. 



Heavy quark states can be an important probe of the dynamics of the Quark Gluon Plasma (QGP). While charmo- 
nium suppression [1] has been observed at a range of energies, a number of sometimes competing effects complicate 
the interpretation of these patterns. Many of these difficulties are not present for bottomonium making suppression pat- 
terns observed in that regime in principle more straightforward to interpret. Intriguing early results from experiments 
at CMS [2] and STAR [3] demonstrating suppression in the Upsilon system suggest a lattice calculation is timely and 
may help to shed light on the link to the spectrum of bound states. 

In this paper we use dynamical anisotropic lattice gauge configurations and treat the bottom quarks using non- 
relativistic QCD (NRQCD) at a range of temperatures, 0.427; <T < 2.097;. NRQCD has been used extensively to 
study heavy quark physics both on and off the lattice. At finite temperature this formalism offers a number of technical 
advantages over the relativistic theory. Firstly, the temperature dependence in NRQCD correlators is due only to the 
thermal medium, with no contribution from thermal boundary conditions. Secondly, there is no nontrivial spectral 
weight at zero energy, which would yield a constant time-independent contribution to the correlators complicating the 
analysis [4] and casting doubt on the results for melting or survival of charmonium at high temperatures [5]. 

The non-relativistic theory is also the first effective theory to be obtained when integrating out the UV degrees of 
freedom [6, 7] and requires only that the heavy quark mass, M^T which is reasonable for ^7-quarks at the temperatures 
in our simulations, up to 27; ^ 400 MeV. A more detailed discussion is in Ref. [8] and references therein where some 
results shown here have also appeared. A complete MEM analysis will be presented in a forthcoming paper. 



In this paper we present results from a study of the temperature dependence of S and P waves (in the T and Xhi 
channels respectively). From Ref. [7] the correlators for S and P states, in continuum NRQCD with Ep = p^ /2M and 
in the absence of interactions, are 
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indicating a power-law decay at large euclidean time. Although in a lattice simulation this behaviour can be modified 
by interactions and lattice artifacts it provides a useful guide for what to expect. 

In this paper we also present preliminary results from a maximum entropy method (MEM) [9] analysis of bottomo- 
nium S waves. The dataset used in the MEM analysis is modified slightly - including more configurations and a larger 
range of temperatures as well as a more accurate nonperturbative tuning of the bare anisotropy, ^ = ajax in the action. 
The gauge configurations are those used in Ref. [11] with ^ = a.v/«T = 6. NRQCD bottomonium (point) propagators, 
using an action including terms up to ff{v^), were computed and analysed for a number of S and P wave channels. 
The analysis presented in [8] and described briefly below was based on a subset of the temperatures now available (as 
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Some lattice details for this work. The lattice spacing 


is a^' fa 7.35 GeV, set using the IP- IS 


spin-averag 


ed splitting in charmonium [12]. 
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500 
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408 
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441 [458] 


2.09 


1000 


500 



indicated in the table). Full details of the nonperturbative tuning of the anisotropy are in Refs. [12, 13]. 



Temperature dependence of the T and Xb states 



To distinguish bound and free states we consider the effective mass, meff(T) = — log q'^}1\ ■) and define a new 
quantity: the effective power, 7eff(T) = _^ G(T+ni)-G(T-nT) ^ p^^. g^^^gg jjj hadronic phase we expect to see the 
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FIGURE 1. The left panel shows the effective mass for the S and P wave states. The right panel is the effective power for the 
same states. Each pane shows a range of temperatures. 

exponential decay characteristic of bound states, with mgff (t) becoming a constant at large t. When quarks are unbound 
we expect the correlators will show a power-law behavior and that 7eff (t) will tend towards a constant. These quantities 
are shown in Fig 1 . We note that the P wave behaviour at T > 71- rules out pure exponential decay, which we interpret 
as a medium modification at temperatures above Tc- In addition, we see a tendency of the P wave to flatten out and 
approach the free continuum behavior at large T (~ 2Tc). 



Maximum entropy analysis. To study in detail the temperature dependence of the bottomonium correlators we 

calculate the spectral functions with MEM using the relation, G(t) = / e^'^'^p (to') , where (O = 2M + co' and 

J-2M 2n 

dropping terms suppressed when M [7]. For the MEM to give reliable results a sufficient number of points in the 
euclidean time direction is required: ^(10) independent lattice points is a good guide. At temperatures close to 2Tc this 
would necessitate a lattice spacing a ~ 0.025 fm, which on an isotropic lattice is prohibitively costly. The anisotropic 
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FIGURE 2. T spectral functions, normalised by the heavy quark mass; at zero temperature (left) and nonzero temperatures (right). 



lattice provides an excellent solution: with ^ a sufficient number of points for MEM to work can be simulated 
for reasonable numerical cost. The MEM analysis is performed using Bryan's algorithm [10] and quadruple numerical 
precision. In the results shown here we have focused on the S waves. We begin by looking at the zero temperature 
spectral functions in the Y (vector) channel. Fig. 2 (left panel) shows the spectral function at the lowest temperature as 
a function of the energy. The dotted lines indicate the ground- and first excited state determined using exponential fits 
to correlators. We observe that there is very good agreement with the MEM spectral function and that the ground- and 
first excited states are clearly determined. The right panel in Fig. 2 shows the temperature dependence in the T channel 
with each pane showing two neighboring temperatures. The ground state appears to survive to the highest temperature 
in our simulations (^ 2Tc), although the height of the peak is diminished. The first excited state is gradually being 
suppressed and is no longer discernible aiT /Tc ^ 1 .68 indicating a melting temperature below 1 .68rc- 



Summary 

In this paper we describe recent results from a study of bottomonium at non-zero temperature from lattice QCD. The 
NRQCD formalism offers a number of advantages at finite temperature making it a promising line of investigation. 
We find evidence of strong temperature dependence in the P waves at T > and indications that while the IS states 
survive up to 2Tc the 2S states are affected much more and appear to have melted at temperatures above 1 ATc- 
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